Interface roughness, valley-orbit coupling and valley manipulation in quantum dots 
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We present a systematic study of interface roughness and its effect on coherent dynamical processes 
in quantum dots. The potential due to a sharp, flat interface lifts the degeneracy of the lowest energy 
valleys and yields a set of valley eigenstates. Interface roughness is characterized by fluctuations 
in the location of the interface and in the magnitude of the potential step. Variations in the 
position of the interface, which are expected to occur on the length scale of the lattice constant, 
reduce the magnitude of the valley-orbit coupling. Variations in the size of the interface potential 
step alter the magnitude of the valley-orbit coupling and induce transitions between different valley 
eigenstates in dynamics involving two (or more) dots. Such transitions can be studied experimentally 
by manipulating the bias between two dots and can be detected by charge sensing. However, if the 
random variable characterizing the position of the interface is correlated over distances of the order 
of a quantum dot, which is unlikely but possible, the phase of the valley-orbit coupling may be 
different in adjacent dots. In this case tunneling between like and opposite valley eigenstates is in 
effect a random variable and cannot be controlled. We suggest a resonant tunneling experiment 
that can identify the matrix elements for tunneling between like and opposite valley eigenstates. 



I. INTRODUCTION 

Of all semiconductors investigated in quantum compu- 
tation (QC) only C, Si, and Ge in group IV have isotopes 
with zero nuclear spin, allowing long coherence times, 
with Si long known as an outstanding candidate.'^'QC ar- 
chitectures are based on donor electron or nuclear spins in 
Si:P,^ ^ single electron spins in gate-defined quantum dots 
(QDs) in Si/SiG«PElor Si/SiOa,*^"^ and SiGe nanowires,^ 
and singlet-triplet qubits.^i^ Experimentally, significant 
progress has been made recently in Si:P donor-based 
deviceJi^l^ and gate-defined Si quantum dots.^ 321 
Graphene has also been proposed for QC,l^and observa- 
tion of the Coulomb blockade effect in a graphene quan- 
tum dot has been reported.'^^ Carbon nanotube quantum 
dots offer another physically interesting possibility, with 
promising experimental advances.!^ The theory of QDs in 
group IV materials has also taken off in recent years.l^^l^ 

The conduction band of group IV semiconductors con- 
sists of a series of equivalent valleys: bulk Si has six, 
Ge has four, whereas in C bulk diamond has six valleys 
while graphite, graphene and nanotubes have two. This 
multiplicity is known t o affec t spin quantum computing 
schemes in Si: P don ors,!^^ Si QDsPSHSni and carbon 
nanotube QDs.'^^'^ Confinement partially lifts the val- 
ley degeneracy, and it is often necessary to consider only 
two valleys. In Si (001) the ^-confinement lowers the en- 
ergy of the z-valleys, perpendicular to the interface, by 
several tens of meV relative to those in the plane of the 
interface. The interface potential produces a valley-orbit 
coupling A between the two £- valleys,^ which is sample- 
dependent .'^Measurements of A have been performed on 
quantum Hall system^^SM] magnetic fields.'^ 

A large valley splitting was recently reported in Si/SiOjaa 



and work is underway to confirm this claim. 

The valley-orbit coupling has been studied theoret- 
ically within th e effect ive mass approximation (EMA) 
by many groups. '^'^^'^IHlll j^gf, |^ determined A using a 
hard-wall boundary, requiring the wave function to van- 
ish at the interface, whereas in Refs. I27I43I44] the wave 
function was allowed to leak into the insulator. Ref. [17] 
highlighted the importance of interface miscuts in reduc- 
ing A. Beyond the EMA, Ref. 45 performed atomistic 
calculations for Si/SiGe, finding that the presence of dif- 
ferent bonds at the interface suppresses A. Refs. 
also included models of interface roughness. Neverthe- 
less, EMA studies to date have not considered systemat- 
ically the effect of interface roughness on the valley-orbit 
coupling and on the energy-level spectrum of multivalley 
semiconductor QDs, highlighting the importance of de- 
vising a systematic way to quantify interface roughness 
and study its interplay with the valley degree of freedom. 

Interface roughness embodies variations in the posi- 
tion of the interface and variations in the magnitude of 
the potential step at the interface. Variations in the po- 
sition of the interface are characterized by two length 
scales. The rms fluctuations are a measure of the verti- 
cal extent of the roughness, while the correlation length 
contains information about its horizontal spread. The 
study of roughness began with the t heore tical work of 
Prange and Nee,!^ followed by Ando,'^SIIil yet theoreti- 
cal and numerical studies have focused mostly on rough- 
ness effects on transport .l^^l^ Experimentally Goodnick 
et aPS] used TEM to study roughness at the Si/Si02 in- 
terface, determining the magnitude of the fluctuations 
and their correlation length. Ourmazd et al^^ also using 
TEM, determined that the transition layer between Si 
and Si02 spans a distance of approximately SA, in other 
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words approximately two monolayers. The rms fluctua- 
tions in the position of the interface are generally found 
to be of the order of 2-4 AfSHm while their correlations 
extend over a few lattice constants. Experimental stud- 
ies indicate that state-of-the-art technology can be used 
to redu ce th e rms roughness in MOSFETs to the range 
0.7-1 A.ISSISI] Later experiments have hinted that the cor- 
relation length characterizing roughness may be larger. 
Yoshinobu et used atomic force microscopy to infer 
a correlation length of 15 nm, a conclusion upheld by 
Gotoh et aP^ using scanning tunneling microscopy, who 
found a correlation length of 23-26 nm. Even small am- 
plitude rms fluctuations can result in short time scales 
for mixing valley-split levels, since the wave vectors of 
the bottoms of the two valleys correspond to very short 
length scales (sa lA in Si). All the above works have con- 
centrated on variations in the position of the interface. 
Variations in the size of the potential step, that is, the 
local roughness, have not been studied to date. 

In this article we extend the EMA to study the ef- 
fect of interface roughness on multivalley QDs. The ap- 
proach adopted in this work is equivalent to that of Ref. 
1311 except that an analytical wave function is used here 
rather than one determined numerically. Our theory is 
applicable provided that fluctuations in the location of 
the interface are much smaller than the spatial extent of 
the 2DEG, and the fluctuations in the magnitude of the 
interface potential step are small compared to its aver- 
age value. We demonstrate that interface roughness has 
two main effects: it suppresses the valley-orbit coupling 
and gives rise to tunneling between valley eigenstates in 
certain dynamical processes in multidot systems. Our 
starting point is the effective mass approximation in the 
presence of valleysj^'^'-with a modified Hamiltonian incor- 
porating interface roughness. Our model of roughness 
accounts for variations in the position of the interface 
and variations in the magnitude of the interface poten- 
tial step. For one electron on one dot [leld) we study 
the way interface roughness enters the expression for A. 
If roughness correlations occur on a spatial scale smaller 
than that of the dot, variations in the position of the 
interface reduce the magnitude of A. For one electron 
in a double quantum dot {le2d), under these circum- 
stances A is expected to be the same on each dot. Local 
roughness can cause A to be different on the two dots 
and induce transitions between valley eigenstates. In- 
terdot transitions between valley eigenstates can be un- 
derstood as coherent rotations of valley states, and offer 
the possibility of manipulating valley eigenstates by ac- 
cessing the appropriate region of the energy spectrum. 
We subsequently study transitions between two-electron 
states in a double quantum dot (2e2d), which represent 
the most promising scenario for observing such transi- 
tions. For roughness correlations of the order of the size 
of the DQD, A is not suppressed, yet its phase may differ 
between adjacent dots. Tunneling between valley eigen- 
states is enabled in this case. The size of this tunneling 
matrix element cannot be controlled, although it can be 



measured experimentally. The findings of this work are 
general, yet when concrete examples are called for we 
will discuss QDs made in a Si/Si02 MOSFET structure, 
which are actively studied in quantum information. 

It should be pointed out that the long spin coherence 
time in Si and the possibility of further enhancing the 
spin coherence time using isotopic purification, which 
eliminates the ^^Si nuclei, give Si an enormous advan- 
tage for quantum information processing. In addition, 
the dominant Si-based semiconductor technology allows, 
in principle, for scalability with the real hope that once 
a few Si spin qubits are demonstrated in the laboratory, 
many qubits could be developed without great difficul- 
ties. These (i.e. long spin coherence time and scalabil- 
ity) are the reasons that Si spin qubits remain a most 
active area of research in the context of quantum com- 
putation architecture studies. However, the conduction 
band valleys in Si, which are degenerate in the simplest 
effective mass approximation, cause considerable difficul- 
ties in fabricating Si spin qubits since the electron spin 
by itself does not define the state of a Si electron. This 
is the context in which our detailed theoretical study in 
this work has been undertaken. The ultimate goal is to 
help understand the role of valleys in Si spin quantum 
computation, in particular, to see if the interface rough- 
ness naturally present at the Si surface could somehow be 
utilized to circumvent or even control the Si spin qubit 
valley degeneracy problem. 

The outline of this article is as follows. In Sec. HIl we 
consider the case of a single quantum dot, valleys and 
valley-orbit coupling. We provide an analytical scheme 
for calculating the valley-orbit coupling. In Sec. |III| we in- 
troduce a simple model of interface roughness and study 
the additional terms it introduces into the leld prob- 
lem. In Sec. |IV| we introduce a model of a double quan- 
tum dot. In Sec. VA we study intervalley transitions 



induced by interface roughness for le2d strong tunneling 
and le2d strong roughness. In Sec. |VB| we study transi- 
tions among 2e2d singlets and among 2e2d triplets, while 



Sec. VI concentrates on the case when interface rough- 
ness correlations are of the order of the size of a quantum 
dot. An extensive discussion of the results, as well as nu- 



merical estimates, are contained in Sec. VII We end with 
a summary of our findings. 



II. SINGLE QUANTUM DOT 

The EMA assumes the overall potential experienced by 
the charge carriers is separable and smooth. Convention- 
ally the EMA is understood to apply to physics parallel 
to the interface, within the 2DEG, as long as the system 
under study is a few lattice constants away from the in- 
terface. The smoothness requirement can be relaxed and 
the EMA generalizecP^ to account for the presence of 
an interface by considering a z-dependent effective mass, 
smoothly interpolating between the two effective masses 
on two sides of the interface plane, as we do in this work. 
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In this section we introduce our theoretical model and 
study a smooth interface as a benchmark for subsequent 
comparison. 



The effective mass for SiGe is the bare electron mass 
mo, w hile f or Si02 it is found to range between 0.3 and 
0.4mo.'^^'^For 'tjj{z) we use a modified Fang-Howard vari- 
ational wave function.'^ With the interface at z = 



A. Model Hamiltonian 

We study first one electron on one dot. The full Hamil- 
tonian Hf)(x,y,z) includes the kinetic and potential en- 
ergy contributions. The confinement potential is ap- 
proximated as quadratic in the xy-plane, while in the 
i-direction one has the interface potential V{z) and elec- 
tric field F 



VD{x,y,z) 



2m* or 



y 
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+ V{z) + eFz. 



(1) 

The dot is located at (x£),0,0) and has a Fock-Darwin 
radius a. The interface potential V will be discussed in 
Sec. |HI[ Within the effective mass approximation (EMA) 
the wave functions for the quantum dot are written as 
D^{x,y,z) = (j)D{x,y) ip{z) U(^{r) e^'''^^ , where, as in our 
previous work, we have introduced a valley index ^ = 
{z, z} that will become relevant below, and fc^ ^ = zL fep = 
±0.85 (27r/asi), with asi the lattice constant of Si.l^aMI 
For the envelopes (j^oix, y) we use the Fock-Darwin states 



(t>D{x,y) 



(2) 



ij^iz) = Me^,z < 



(4) 



N {z + Zo)e^^,z > 0. 



We will also use b = 1/ksi as a length scale characteris- 
tic of the size of the 2DEG along the growth direction. 
Continuity of the wave function at the interface requires 
M = Nzq, and that of (l/rriz) dip/dz gives 



N 



2zo 



(5) 



The parameter zq is given by 



(6) 



The effective mass equation for motion in the z- 
direction, perpendicular to the interface, is 



2m- 



dz^ 



+ V(z) -I- eFz 



^{z)=eMz). (3) 



In this work ki, is fixed at ki, = y ™2 , leaving as 

the only variational parameter. The energy that is to 
be minimized is given by 



e. = {2b^ + 2hzo - z^) + (66' + 4bzo + z^) 



8mi,b kb 



2kl 



r 



(7) 



which includes a factor of 1/2 in the electric field term 
to account for double counting of electrons.!^ 



B. Valley-orbit coupling for a sharp flat interface 

The conduction band of some materials consists of a 
series of equivalent valleys. In such circumstances, for 
a smooth potential the effective mass equation has two 
equivalent solutions. If the potential is not smooth the 
solutions could be mixed, leading to valley-orbit coupling. 
Consider the basis states {D^}, where D — L,R repre- 
sents either of the left and right dots. The z and z states 
have a vanishingly small overlap, which is neglected here. 
With this approximation these states are orthogonal. To 



determine the one-dot energy levels we require first the 
matrix elements of the Hamiltonian in the basis spanned 
by {DJ = either {LJ or 



Id 



£d 







A 



D 



Az5 





(8) 



where ££> is the confinement energy, which is large, thus 
we consider only the lowest level = generally of 
the order of tens of meV. The lowest energy level for each 
dot consists of two valleys connected by a valley-orbit 
coupling A^) — |A£)|e~*'^". The eigenstates of ifieid 
can be expressed as 



1 

7^ 



\D,)±ef'"\D,)). 



(9) 
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In this work intervalley refers to matrix elements con- 
necting the z, z states, whereas between valley eigenstates 
refers to matrix elements connecting + and — states. The 
valley eigenstates are the same for both dots if the inter- 
face is sharp along the growth direction and flat perpen- 
dicular to it. The valley-orbit coupling is represented by 
the matrix elemenli^ 

^D = {D^ly + eFz)\D,). (10) 



Its magnitude gives the valley splitting lA/jj. The dom- 
inant contribution to comes from the interface po- 
tential, while the contribution due to F has been found 
to be negligible.'^ For a perfectly smooth and perfectly 
sharp interface A^i is equal to 



Ad= Uo f f dxdy\cl){x,y)f f dz Q{-z) \ij{zf e~^'''^' uHr) u,{r) ^ Uo z^^ ^ 



(11) 



where ^ ki, + iQ^ — 2ik^. The matrix element between the same valley states = {D^\V\D^) is given by 



Ad = Uo J Jdxdy\<j>ix,y)f JdzQ{-z)mz)f\u^{r)\^ = UoN^z^oT. 



K,Q^ 



Pi 



(12) 



where = fcb + iQz ■ We have used the general nota- 
tion ^ for the valley indices in A^j since we anticipate 
{D^\V\D^) = {Di\V\Ds). Detailed derivations of Eqs. 
(Ill and (12) are presented in Appendix [A| The coefh- 



cients c]^ for Si were determined in Ref. [55] Using these 
values, we find that by far the biggest contributions to 
A/3 and A^j come from the terms with Qz = 0. Given 
that the magnitude of \kb—2iko\ « 2fco, it follows that A^ 
can be approximated simply by A/j « iUo S A^^ 2:g/(2fco), 
where S = c^'^jc^- The Umklapp terms, which have 
Qz = 2^ (27r/asi), give imaginary contributions and sum 
to zero. We will refer to the valley-orbit coupling aris- 
ing from the sharp flat interface as the global valley-orbit 
coupling Aq, without the dot subscript, and similarly for 
the global A^i = Aq. 

To obtain numerical estimates of |Ao| in Si/Si02 and 
Si/SiGe an interfacial electric field of 150 kV/cm is as- 
sumed, which is the same value that was used in Ref. 
1311 in order to enable a direct comparison. For a Si/Si02 
interface, with Uo ~ 3eV (as discussed in Ref. [51]) and 
mjj = 0.4mo, b is optimized at 1.058 nm, for which we 
find |Ao| «0.11 meV. Since the effective mass is not well 
deflned for Si02, we have also used to;, — O.Stoq, find- 
ing b « l.Glnm and |Ao| «0.G8 meV. For a Si/SiGe 
interface, with Uo ~ 150meV, b is optimized at 1.215 
nm, whence |Ao| ~0.1 meV. Considering the relatively 
simple variational function that is used in this work, 
the comparison with accurate numeric^^ is extremely 
encouraging. The real and imaginary parts of Aq are 
given to a very good approximation by kh / \/k1 + 4fcg 
and 2kol yjkl + Ak'^. Consequently, with ko fixed, as the 
interface potential step becomes stronger the imaginary 
part of Aq also increases. For the Si/Si02 interface the 
ratio of the imaginary and real parts of Aq is found to 



be 3.5, while for the Si/SiGe interface it is 15.7.E21 The 
ratio of the off-diagonal matrix element Aq of the inter- 
face potential to its diagonal matrix element Aq is given 
approximately by kh/{kh — 2iko), and the magnitude of 
this ratio is 0.28 for the Si/Si02 interface and 0.06 for the 
Si/SiGe interface. The final result varies significantly as 
a function of b^ 



III. INTERFACE ROUGHNESS 

An interface may be observed to be sharp and of a high 
quality, yet the transition layer between two materials 
still always has a finite width. The atoms on the two sides 
of the interface have different bonding, as do the atoms 
at the interface. For example between Si on one side and 
Si02 on the other there exists a region a few monolay- 
ers thick composed of SiO^;, with < a; < 2.^^ Conse- 
quently, even when the interface is sharp the finite tran- 
sition layer alters the energy landscape.!^ Experimen- 
tally rough ness c an be detected by transmission electron 
microscopy,'^Sl52l atomic force microscopy,'^ scanning tun- 
neling microscopy,'^ and other surface characterization 
techniques such as Auger spectroscopy, angle-resolved 
photoemission spectroscopy, low-energy electron diffrac- 
tion and resonant high-energy electron diffraction!'®^ 

The notion of a sharp barrier breaks down at the 
atomic scale. The ideal method for dealing with this 
issue is to use the exact microscopic structure of the in- 
terface. Such a treatment is only possible numerically, 
yet even then the interface structure is not known exactly 
and, when one of the materials is amorphous, it can only 
be dealt with approximately and/or statistically. There- 
fore, there exists a scope for gaining analytical insight 
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by refining the effective mass approximation to include 
the effects of interface roughness. Within the EMA the 
interface is modeled in such a way that the atomic struc- 
ture is smoothed out, and an effective mass appropriate 
for the insulator is used. In the same way as the smooth- 
ness requirement can be relaxed as long as the region 
under study is a few monolayers away from the interface, 
one can also relax the requirement that the interface be 
sharp and flat. Corrections to this picture can be made 
so that the location of the interface varies on an atomic 
scale but with an amplitude much smaller than the ef- 
fective size of the confinement. Then the xy-motion and 
the z-motion are still approximately separable. Beyond 
this refinement, a random potential is added to allow the 
magnitude of the interface step to fluctuate. These two 
features can be regarded as corrections to the effective 
mass picture, which is entirely based on bulk parame- 
ters. 

A si mpl e model, which was i ntrodu ced by Prange 
and Ne^^ and developed by Ando,'^2lII] has become the 
standard model used in studies of surface roughness. 
Within this model the interface is described by a po- 
tential V[z ~ (^(x, y)], in which the fluctuating position of 
the interface C(a;, y) is treated as a random function of x 
and y. Prange and Ne^^ assumed that the probability 
for an electron to be specularly reflected at the interface 
is nearly unity, with a boundary condit ion th at the wave 
function vanish at the interface. Andd^^*^ approached 
the problem by expanding V in C and retaining the term 
linear in This approach may be satisfactory for cal- 
culating the scattering matrix element due to roughness. 
Yet, as Refs. ISOlSTI themselves clearly state, there is no 
rigorous justification for this expansion. As we will show 
in this work, indeed this expansion is not at all justified 
when the effect of roughness on the valley-orbit coupling 
is required. The breakdown of the perturbative expan- 
sion is most emphatically manifest in Section |VI[ when 
long interface correlations are considered, yet the expan- 
sion is also not valid if the ^-correlations span a small 
spatial scale, as Section [V] demonstrates. 



A. Model of roughness 

We treat interface roughness using a slightly different 
approach from that of Refs. 50 and 51, We consider an 
interface potential given by 

Vix,y,z)^ [Uo + V{x,y)]f[z-ax,y)]. (13) 

where Uq ^ V ioi all x and y. We take a general func- 
tion f{z — to describe the interface, where f{z) is a 
well-defined function of z, while C is a random function 
of X and y accounting for the variation in the location of 
the interface. For an interface that is perfectly sharp in 
the i-direction and perfectly smooth in the xy-plane (the 
ideal interface), one can write V(z) = Uq Q{—z). For an 
interface whose location fluctuates in the xy plane, but is 



still sharp, one can use f{z — C) = ©[— (-^^ — C)]- The po- 
tential V{x, y) represents the local roughness, which gives 
fluctuations in the magnitude of the potential step at the 
interface. We treat V{x, y) perturbatively, since we ex- 
pect V{x, y) ^ C/qi but we do not do perturbation theory 
in C,{x,y). In this problem the criterion for smallness is 
generally nontrivial. The discussion of this criterion will 
follow naturally after the analytical expressions for the 
effect of the roughness are presented. 

Aside from variations in the location of the interface, 
sources of roughness include dangling bonds, impurities, 
and variations in bond length in the vicinity of the in- 
terface. All of these contribute to the local roughness 
V{x,y). We expect, however, that variations in bond 
length will occur on the atomic length scale, which is or- 
ders of magnitude smaller than the spatial extent of a 
quantum dot. This implies that the effect of bond length 
fluctuations averages to zero over the quantum dot and 
is washed out. We do not consider it explicitly in what 
follows. 



B. Applicability of the EMA 

The location of the interface C is a random function 
of x and y. It has an average value {(,) — 0, while its 
variance will be denoted simply by (C^)- The notation 
denotes the average over the xy-plane, while ()qd de- 
notes the average over x and y spanning the region oc- 

cupied by the quantum dot = l/(7ra^) / / dxdye . 
It is assumed that the fluctuations around {Q are small 
compared to 6, which represents the spatial scale over 
which the wave function changes by a significant amount 
in the z-direction (i.e. the size of the confinement.) The 
spatial scale of the roughness divided by the spatial scale 
over which the wave function changes significantly is an 
important parameter. The size of this parameter in com- 
parison to unity determines the applicability of the EMA. 
It is necessary for C,{x, y) to be small, so that the poten- 
tial acting on the charge carriers remains approximately 
separable and use of the effective mass approximation is 
justified. If 1(^1 <C a, 5 the fluctuation in the z-direction is 
considerably less than the spatial extent in all three di- 
rections, so that the wave function is still approximately 
separable in x,y,z. This condition must be augmented 
by a formal criterion according to which the potential is 
weak. In principle all applied electric fields, such as F in 
this problem, violate the criterion for being weak, since 
the potential eFz eventually becomes large on any scale 
relevant to the problem once z exceeds a certain thresh- 
old. This fact however is not usually a problem as long as 
the potential V does not vary too fast on lattice length 
scales, that is, as long as eFasi is significantly smaller 
than the average atomic potential in a unit cell. To sum- 
marize, the EMA should remain valid as long as the ran- 
dom potential V{x,y) ^ Uq for all x,y and the fluctu- 
ations in the location of the interface are much smaller 
than the spatial extent of the 2DEG. 
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C. Roughness effect on the valley-orbit coupling 

We use the notation Ao for the unperturbed valley- 
orbit coupling in the absence of interface roughness, 
which is the same on both dots. The valley-orbit cou- 



pling in the presence of roughness is denoted by A. We 
require |A — Ao| ^ |Ao| for any x,y. In the presence of 
an interface located aX z — C, the wave function must be 
altered as ~^ — ()■ The matrix elements of the 
interface potential are found in the same way as before 



An = 



dxdy(C/o + y)|0(x,y)|2 / e[-(z - C)] |^(^ - C)l' 



u*(r)u_^(r) = 



V) 



QD 



(14) 



The effect of the interface is contained in two terms: the 
spatial average of the exponential e*"? , and the spatial 
average of the product y) e"'^ ^. We recall that in 
Si the contribution to A due to one wave vector, namely 
Qz = 0, overwhelms all the others. Therefore a good 
indication of the effect of interface roughness is given by 
the averages (e^*''"^)Q_D and {V{x,y) e^*'"' '^)q£i- The ap- 
pearance of the simple exponential oscillatory terms such 



as e"*5 ^ is a consequence of the fact that we are using the 
effective mass approximation. In a more complicated nu- 
merical treatment we still expect a strong dependence 
on ^, though it will likely not be expressible as a sim- 
ple exponential. Since it is predominantly the term with 
Qz — that contributes to A/^, this quantity should not 
be affected by the presence of C- This can be seen imme- 
diately by considering 



Ad= / dxdy{Un + V)\(j){x,y)\^ dzQ[~{z 



We distinguish between V{x,y), which can be treated 
perturbatively, and e*''^^, which is non-perturbative. 
The most suitable approach to the problem is determined 
by the spatial scales involved. Experiments indicate that 
the spatial scales characterizing ( are 2 — AA for the rms 
and 5 — ISA for the correlation length. Several studies 
have indicated that the roughness in Si/Si02 could be 
reduced, to an rms of the order of « 1 A^^ and even as 
far as O.lA^ At the same time, a series of experiments 
using scanning tunneling microscopy suggested that the 
horizontal correlations of C could stretch over distances 
of the order of 30nm. We will use short correlations to 
refer to spatial correlations of C on a scale much smaller 
than the typical extent of a quantum dot (approximately 
50nm), and long correlations for spatial correlations of 
C on a scale of the order of or greater than the typi- 
cal extent of a quantum dot. The problem of valley-orbit 
coupling is vastly different for short interface correlations 
and long interface correlations. Physically, fluctuations 
associated with a random interface ought to occur on the 
scale of the lattice constant. We anticipate that the true 
physical picture is represented by short correlations, but 
given the existing uncertainty we will study both cases 
in this paper. 

In the most general case the valley-orbit coupling may 
be different on the two dots. Since the valley-orbit cou- 
pling is a complex number, in principle one should expect 



0] - OP \u,ir)\' = ^'^'o{iUo + V))Qn _ ^^^^ 

I 

both the amplitude and the phase to be different between 
the two dots. This fact is not of much significance when 
a sharp flat interface is considered, yet is of enormous im- 
portance when interface roughness is concerned. In ad- 
dition, roughness is expected to show significant sample 
dependence. In double quantum dots in different samples 
the roughness profile on the left dot may be qualitatively 
similar to or very different from the roughness profile on 
the right dot. One also expects a significant difference 
between Si02 and SiGe, with the former generally hav- 
ing larger interface roughness. From this perspective it 
will be useful to be able to compare e.g. gate-defined 
quantum dots with disordered quantum dots. 

In addition, we have noted that, in sums of the type 
J2c cce*^ '', the only sizable terms correspond to the re- 
ciprocal lattice vector G = (000). We infer that terms 
with wave vectors larger than that, corresponding to spa- 
tial scales of lA or less (atomic scale roughness) , should 
not be important in determining Ajj and A^i. This rein- 
forces the observation of Sec. |III A| concerning variations 
in bond length. 



IV. DOUBLE QUANTUM DOT 

In a single quantum dot interface roughness affects the 
valley-orbit coupling. In a double quantum dot it affects 
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in addition the relative values of the dot-averaged valley- 
orbit splitting Ad on the two dots, as well as interdot 
tunneling, potentially giving rise to interesting electron 
dynamics. The confinement of a double quantum dot is 
modeled by a biquadratic potential 



^r I X ^0 r,.. Ax 



a'^ 



-eEx + eFz + V(z). 

(16) 

As in Ref. [301 we orthogonalize the wave functions 
{D^} to get the Wannier functions {D^}, given by 
Lf = g-R; ^ ^j^^ _ R( gLi ^ where q 



(1 - Vl - and ; = (^el^c)- This ensures that 
{R(\L^) = 0. One diagonalizes the Hamiltonians Ho 
and determines the eigenf unctions {D±]. We will con- 
sider one electron initialized in a DQD {le2d) as in a QD 
charge qubit, as well as two electrons {2e2d), which can 
form singlet and triplet states (given in Ref. [301) Hence- 
forth we work in this basis, and to indicate this fact we 
augment all quantities with a tilde, thus A etc. 

Our principal assumption in setting up the problem is 
that the physics specific to the valleys can be considered 
separately for the left and right dots. The starting point 
of our approach envisages two dots and one electron pri- 
marily localized on one of the dots. The valley degree of 
freedom implies that there are two possible single-particle 
states for each dot, while the interface potential gives a 
coupling between the two valleys on each dot. We diago- 
nalize the single-particle Hamiltonian separately for the 
left and right dots, and obtain the valley splitting and 
valley eigenstates for each dot. We neglect matrix ele- 
ments of the interface potential connecting the left and 
right dots, which represent terms of order I A. Under 
these circumstances (/>_d is not a function of position, a 
fact that provides a crucial simplifying assumption in our 
algebra. 

The assumption can be recast into the statement that 
there is no interdot tunneling between different valleys, 
meaning that the tunneling matrix element only connects 
states from the same valley. One can tunnel between \L^) 
and \R)^, but not \L^) and 1^?)-^. Tunneling between 
valley eigenstates is determined only by the phases of 
each valley wave function in the valley eigenstates on each 
dot, which in turn are determined only by the details of 
the interface potential in the region of the dot. 

A full treatment of the le 2d and 2e2d problems would 
take as its starting point the four valley wave functions 
{|Lj)} and It takes into account all the ma- 

trix elements of the interface potential, including ma- 



trix elements of the form (L^|V|i?_^). In this case the 
eigenstates are mixtures of all four functions and 
a far longer and more cumbersome approach. 
Since for a typical DQD Z <C 1 by several orders of magni- 
tude, terms of order I A are unlikely to play an important 
role, and it is reasonable to anticipate that our simplified 
approach contains the dominant physics. Given that A 
differs from A by terms of order ^A, which are neglected 
here, we can consider A w A. 



A. One-electron Hamiltonian 

The most general Hamiltonian for one electron in a 
DQD in the basis {!/+, R^, i?-} is 



Hle2d — £^ 



f\K\ 

^0+ + 



-JA^I 



to — h 
-S + \A^\ 



to+- \ 
to— 

-S-\A^\/ 
(17) 
and 



where the detuning is defined as 6 

|A^| = |A| ± v±±- The tunneling parameter tg^ is 

defined as ion ~ {L+\H\R^) . We will use the notation 

-D^, — {I)^\V{x,y)\D^i), with corresponding notation for 
the ± states. Throughout the remainder of this article, 
we will also use v to refer (generally) to the magnitude 
of the off-diagonal term an important parameter 

in determining several regimes of qualitatively different 
physics. 



B. Two-electron Hamiltonian 

For 2e2d, the interdot tunneling parameter t has a 
single-particle par t to an d a Coulomb-interaction (en- 
hancement) part g p8 | 29 | 70 |^]-^g total being i — io + s, with 



S = {L^i^L^i^\Ve,\L]^'R 



(l)n(2)\ 



The Coulomb exchange term j is much smaller than 
the other terms appearing in the Hamiltonian and is 
not needed at all in this paper. We also define the 
critical detuning as Sc = {u — k), where u repre- 
sents the on-site Coulomb interaction and k the di- 
rect Coul omb interaction between electrons on the two 
dots.'22122l The effective detuning is defined as S = S — Sc, 
measured from this critical detuning. In the basis 
{S^R^S^R^S^R^S^R^S^R^S^l^^S^R} used in Ref. M 

the most general singlet Hamiltonian is Hs — Hg'^ + Vg, 
where 
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(1) 



/-2|A| i++y/2 

i++%/2-(5-2|A| r+\/2 

2|A| t++V2 L+ 

t++^/2-~5 + 2\K\i\_^/2 

f+_\/2 i++ 

t_+\/2 f++ 

V r+ t++ i++ -I) 



(18) 



i/g^^ was given in Ref. 30 without the tunnehng matrix element Ys represents the matrix elements of the local 

roughness between singlet wave functions, which are discussed in Appendix [B] It is given by 
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^/2^-«^ 











(19) 



The most general triplet Hamihonian in the basis {T^^,Tff,Tf^,T^^,Tf^} used in Ref. iSOjis TJr = Hj, 
where 



(1) 



H. 



(1) 





-2|A| 


















2|A| 








t-+ 
















t++ 
















-t++ 


V 






i++ - 




-s J 



Vt, 



(20) 



and 



Vt = 

























+ + 







-'-^ 




















(21) 



The matrix elements of Vt are discussed in Appendix |B] If one assumes that the valley-orbit coupling has the 

same form in both quantum dots then both have the 
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(C^) is the variance of (. In this case, since 



0.5 1 1.5 

I, rms (Angstrom) 



Figure 1; Valley spliting |A| as a function of the rms fluc- 
tuations in the location of the interface ("rms = \/ (C^)- It is 
assumed that ( is described by a Gaussian correlation func- 
tion that varies on spatial scales considerably smaller than the 
extent of the quantum dot. An optimal value of A of O.lmeV 
has been assumed, as determined in Sec. |IIB| for a sharp flat 
interface. 



same valley splitting and valley eigenstates. When elec- 
trons tunnel between the dots they tunnel between like 
valley eigenstates, and interdot intervalley tunneling is 

strongly suppressed. In other words = above. 

In the absence of any interface roughness, with = 

V = 0, the spectrum of two electrons in a DQD can be 
clearly divided into three branches with different valley 

composition,!^ which may be labeled by ++, ,H — . 

These branches may be mixed by interface roughness. 
As we will discuss in the sections below, two features of 
this mixing are noticeable. Firstly, roughness only mixes 
certain branches of the spectrum. There is no mixing 

between the ++ and branches, but both ++ and 

are mixed by roughness with H — . Moreover, there 

is no mixing of states within any of the + + / /H — 

subspaces by the local roughness. Secondly, in situations 
in which tj 



and V couple 



is nonzero (see Sec. VI 

different states: no states are coupled by both v and 

for one or two electrons in a DQD. Therefore, quite gen- 
erally, variations in the magnitude of the potential step 
and variations in the location of the interface have qual- 
itatively different effects on interdot electron dynamics. 



V. SHORT INTERFACE CORRELATIONS 

If the roughness varies over a length scale of a few A, 
then the average of e^**^" over a quantum dot is equiva- 
lent to an average over the xy-plane, {()qd = (C); which 
in turn is equivalent to an average over the realizations 
of C, that is, it is the characteristic function of ^. If one 
assumes C to be described by a Gaussian function, then 
the characteristic function of ^ is given by e"" ' , where 



dCe 



(22) 



which implies that A is reduced by a factor of e~ 
This is plotted in Fig. [T] The reduction of A is in qualita- 
tive agreement with the conclusions of Ref. 45, although 
the spirit of Ref. |3S| is different from this work and a 
direct comparison is not possible. Similarly, a smeared 
interface potential can mimic roughness, allowing the po- 
sition of the interface to be spread out over a range of 
spatial locations,^'* and resulting in a decrease in A. 

The variance of ({x,y) can be extracted from exper- 
iment, and its magnitude is critical. For y/JC^ = 2A, 
the valley splitting A is suppressed by a factor of e^^ = 
3 X 10^^. On the other hand, for roughness of the order of 
Crms — 0-7 A, the exponential factor is e~°-^^ « 0.38. The 
optimal value of A « O.lmeV that we have found above 
for a sharp flat interface would be reduced to ~ 40fieV, 
which should be visible in a dilution refrigerator. At the 
same time we recall that the numerical calculations of 
Ref. l44l found an optimal A of ~ 0.25meV, which would 
be reduced to O.lmeV when Crms — 0-7 A is taken into 
account. We will use this value in our estimates below. 
Thus far the magnitude and phase of A are the same 

on each dot, so = everywhere in this section 

while to-i-i- — to = ^0 smd t^^ = t = t. The 

term {V{x, y) e^* '^^ ^)qd in turn must be taken as a phe- 
nomenological parameter to be measured experimentally. 
Noting that wf^ 



v^g we find 



5£±(l/2) (^"'-vg + e 



'i<t>L 



(1/2) {e-"l'"vg--"l'°^~° 



(23) 



If A is approximately imaginary so that 
expressions simplify to 



i1± 



7r/2 these 



(24) 



(i/2)Rew 



7.D 

ZZ ' 



It is seen that and therefore Vzzi gives the energy 

offset between two dots that otherwise have the same val- 
ley splitting. Knowledge of which can be obtained 
from experiment, allows this diagonal term, important in 
inducing transitions between valley eigenstates, to be es- 
timated. To summarize, for short interface correlations C, 
gives us the overall magnitude of A, while v gives us the 
difference in A between dots. Below we consider the case 
when the local roughness matrix elements v are much 
smaller than A. 



A. One-electron transitions 

We study first the case of one electron on two dots. 
The physics depends on the relative magnitude of v and 
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to- For to ^ V, we diagonalize the Hamiltonian Hie2d 
with aU the roughness terms set to zero, and denote the 
eigenenergies by A^, A^. In this scenario, the lowest 
energy state, which has energy Af , which is the state 
most easily initialized, does not cross any other state. 
Its closest approach to another state is of the order of 
« to- Since ^ v by assumption, valley manipulation 
by means of interface roughness is not feasible in this 
setup. 

For V ^ to, the bare Hamiltonian is now i?ie2d with 
the local roughness terms but without interdot tunneling, 
in other words to initially set to zero. The eigenstates of 
the bare Hamiltonian are localized on the left or the right 
dot. The eigenenergies are denoted A^, A^. Since to is 
set to zero in the bare Hamiltonian, it is trivial to show 
that the eigenstates of i?ie2d in this case do not depend 
on S. It is assumed that one initializes the state with 
energy A;^ and sweeps the detuning until one reaches the 
point where this level crosses the state with energy A^. 
At this point the tunneling matrix element between these 
two levels is 



2LR _ 



\vl\' 



Af 



(25) 



\vr\' 



The critical crossing occurs at the point where A^, which 
depends linearly on the detuning, equals A;^ , which does 
not depend on the detuning. This condition yields a value 
of the detuning of 



(A^ + A^)2 +4|5f_| 



(A« + A:«)2 + 4|{)f _ 



(26) 

The detuning is kept at this value for an amount of time 
« h/i^^. The state of the electron after this time can be 
probed by sweeping the detuning back to 5 3> 0, then do- 
ing charge sensing using a quantum point contact. This 
will reveal whether the electron has tunneled onto the 
left dot or remained on the right dot. 

For weak local roughness the situation described here 
is simply ordinary single-particle tunneling including a 
roughness correction to the energy eigenvalues. The 
eigenstates depend on v, which is not known a priori. 
Yet the time scale governing the state mixing is given by 
the single- particle tunneling matrix element to, which is 
generally small. The size of this matrix element can be 
further reduced by the barrier gate, thus enabling exper- 
iment to increase the time scale for operations to a value 
which can be reliably handled in the laboratory. 

The above demonstrates that, although intervalley 
transitions are extremely slow on one dot, they can be in- 
duced by interface roughness in a coherent manner. Nev- 
ertheless, operations using the single-electron two-dot 
spectrum suffer from the ambiguous initialization prob- 
lems that hamper attempts to manipulate one-electron 



states. Initialization of -D_ is ambiguous since A^i may 
not be known. Moreover, one disadvantage of to being 
small is that it is difficult to know the relative magnitude 
of to and v a priori. 



B. Two-electron transitions 

For strong tunneling we first diagonalize the Hamil- 
tonian with all the roughness parameters set to zero, 
v±± — w-tzp — i±zp — 0. We obtain in this way the 
basis of eigenstates for the strong tunneling case, which 
coincides with the singlet and triplet states used in Ref. 
1301 We work in this basis and add the roughness terms 
v±± and uj-zp perturbatively. Since interface roughness 
does not mix singlets with triplets, we will study the sin- 
glets and the triplets separately. Since the valley-orbit 
coupling is the same on both dots, we still have = 0. 
All relevant quantities are defined in Appendix [C| 

For strong roughness one would have to diagonalize the 
Hamiltonian with all tunneling parameters set to zero, 
obtaining the basis of strong roughness eigenstates. One 
would then add the intravalley tunneling terms t±± as 
a perturbation. Theoretically, the possibility also exists 
that on one dot v ^ i, while on the other v Yet we 
have noted that the tunneling matrix element t contains 
a contribution due to the Coulomb interaction, denoted 
by s, which greatly enhances its magnitude with respect 
to to- One can obtain insight from recent noise measure- 
ments on Si QDs,'^ in which dangling bonds give rise 
to fluctuations in the confinement and barrier potentials. 
Such noise measurements give an indication of the mag- 
nitude of the potential fluctuations one can expect from 
dangling bonds, as well as random impurities, which both 
contribute to v. The size of the fluctuations has been 
measured to be « 0.45^eV^ whereas we expect i to be 
of the order of tens of /leV^ The case v ^ i is therefore 
much less likely in this setup and we do not consider it 
in what follows. 

Since roughness is spin-independent and does not mix 
singlets with triplets, we discuss the singlet and triplet 
states separately. The bare singlet and triplet energy 
levels in the absence of roughness have been calculated 
in Ref. 1291 For convenience we have re-plotted them in 
this paper, showing the singlet and triplet levels sepa- 
rately in Figs. [2] and [3j To determine which levels can 
be mixed by roughness we need to identify pairs of levels 
that cross and find out whether a nonzero matrix ele- 
ment of interface roughness can couple them. The sin- 
glet energy levels for i ^ u are plotted in Fig. |2] Here 
it is more convenient to use the dimensionless detuning 
as 6/{2deo), where d = Xo/a is the dimensionless half- 
interdot-distance. As Fig.[2]shows, the lowest energy sin- 
glet levels do not cross. Even in the case of large A, when 
the lowest singlet state is easily initialized, there is no 
convenient point in the singlet energy level spectrum at 
which one can controUably induce mixing between valley 
eigenstates with different valley composition. The singlet 
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Figure 2: Singlet levels in a double quantum dot with a valley 
splitting of O.lmeV in the absence of interface roughness, as 
determined in Ref. [29l 



0.6 
0.4 



> 0.2 
o 

>, 



iS -0.2 



-0.4 
-0.6 



-■^^ R,R- 






" tfR* 




UR. 


L,R.+/-L.R, 




- - _ L,R.+/-L.R, 


- L R. 




LR 






"-^^R.R. 


1 


, , , 1 


, , , 1 , , , , 1 



0.55 0.56 0.57 0.58 

Dimensionless detuning 



0.59 



Figure 3: Triplet levels in a double quantum dot with a valley 
splitting of O.lmeV in the absence of interface roughness, as 
determined in Ref. 1291 



states for i 3> w are not suitable for valley manipulation 
by means of interface roughness. 

We investigate next the triplet states for t -D. It is 
easiest to access the polarized triplets by starting out 
on one dot (for example R) and applying a large magnetic 
field to lower the energy of the r+ states. This results 
in unambiguous initialization, regardless of the magni- 
tude of the valley splitting, and is easily accomplished 
experimentally since g = 2 in Si. Under these circum- 
stances one needs to consider only the manifold of r+ 
states, which comprises five wave functions. The triplet 
energies are plotted Fig.JS] Far away from the anticross- 
ing (at large detuning S ^ 0) the eigenstates are still 
{f^R^fl^,f^^,f^f,ff'^}. As the detuning is swept 
towards 5^0, the eigenstates evolve into combinations 
of these wave functions, which are given in the appendix. 



Out of the five levels of interest to us, T^:!} evolves into 
r^_, while T:^^ remains unchanged. These two energy 
levels, Tj5_ and T^^ cross at a given value of the detun- 
ing. The critical crossing occurs at the point when the 
energy of , which depends on 6, matches that of T^^, 
which does not depend on 6. Using the energies given in 
the appendix we obtain for the detuning at this crossing 
point S = (2A^ — t'^)/ A. The matrix element mixing the 
two states is 



{T^'1'\V\T<_) = 



t 



.<2 



2^2 



(27) 



Therefore, if one initializes and drives the system 

to this crossing point, the two levels can be mixed by 
interface roughness. The triplet T:^^ is mixed by a small 
amount with T:^^ and Ir+ , but this amount is negligible, 
since it is determined by the ratio \t/S\ <C 1. Therefore 
at the point where the two lowest triplet energy levels 
cross these levels are to a very good approximation 2^;^ 
and T^^. The resulting configuration can be probed by 
charge sensing, since T.^_ goes back to (0,2) whereas 
remains in the (1,1) configuration. This process repre- 
sents a coherent rotation of valley eigenstates. 

We note that for rapid adiabatic passage (i.e. fast 
compared to h/ {T:^^\V\T^_)) the effect of roughness is 
analogous to a renormalization of the spectrum, which re- 
mains qualitatively the same, and valley eigenstates are 
not mixed. However for slow adiabatic passage valley 
eigenstates can be mixed. In effect one can use interface 
roughness to manipulate valleys in the same way the in- 
homogeneous nuclear field was used to manipulate spin 
singlet and triplet states.'^ 



VI. LONG INTERFACE CORRELATIONS 

This case of long interface correlations comprises vari- 
ations in on spatial scales comparable to that of the 
quantum dot. Rece nt ex periments have indicated that 
this may be possible.'^^'^ In order to understand the im- 
plications of this possibility some preliminary considera- 
tions must be brought forth. Primarily among these, it is 
necessary to recall that, if both dots have the same valley- 
orbit coupling, they have the same valley splitting and 
the same valley eigenstates. Under such circumstances 
when electrons tunnel between the dots they tunnel be- 
tween like valley eigenstates. Since this valley-orbit cou- 
pling can vary due to interface roughness, we should in 
principle expect both the amplitude and the phase of 
the valley-orbit coupling to be different between the two 
dots. Since the tunneling matrix element to is defined as 
to - {L\Ho\R), writing Ad = lAole^"^" gives io±± = 
to[l + e*(*-^-'^«)] and io+- = io[l - e'^^^~^'''>]. Since A^, 
does not depe nd exp licitly on position, for 2e2d, using t 
defined in Sec. |iVB[ one has that i±± = i[l + e^('^^--'^«)] 
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and f+_ = i[l - e*('^^-*«)]. The matrix elements of the 
local roughness are unaltered by the relative location of 
the two dots. In all situations discussed in this paper, in 
the case v ^ t, whether the tunneling is intervalley or 
intravalley is a secondary consideration, since tunneling 
is still a perturbation compared to v. The key physics 
is provided by the matrix elements of v, which are the 
same regardless of whether the tunneling is intervalley or 
intravalley. Therefore the peculiar physics of the case in 
which interdot intervalley tunneling is comparable to in- 
terdot intravalley tunneling is manifest when t ^ v, and 
we focus on this regime for the remainder of this section. 
We will concentrate on the more enlightening 2e2d prob- 
lem, in which the most striking consequences of interdot 
intervalley tunneling can be observed. In this section, 

since ^ 0, in the condition t ^ -D it is understood 

that t represents the greater of t and 

There are three qualitatively different situations to 
consider. In the easiest case, correlations stretch beyond 
the size of the entire DQD. The situation is almost iden- 
tical to that of short interface correlations, except the 
value of A is not reduced from its ideal value. If corre- 
lations occur on a spatial scale approximately equal to 
that of one dot, the magnitude of A is essentially the 
same on the two dots but its phase is different: (pL ^ (t>R, 
therefore {(pL — 4>r) ^ 0. From the general definitions 

of t±± and given above, it follows that in general 

both t and are nonzero, and their ratio will de- 
pend on the relative positions of the two dots. Since the 
spatial scale set by 2/co is approximately 0.5A, this ratio 
can vary drastically from one sample to another, and in 
any one sample it is in effect a random number, which 
only experiment can determine reliably. This case cannot 
be treated perturbatively under any circumstances. We 

single out two important possibilities: (a) t = and 

(b) t — t-i In case (a) all the interdot tunneling is 

intervalley tunneling. The situation is qualitatively com- 
pletely different from the case when all the tunneling is 
intravalley, as Figs. |4]and[5]illustrate. None of the singlet 
levels cross other singlets, and none of the triplet levels 
cross other triplets. One cannot accomplish valley ma- 
nipulation with either the singlets or the triplets. Case 
(b) is similar to (a) for the purposes of our paper. As long 
as the intervalley tunneling parameter is appreciable, the 
bottom two triplets are split. In this intermediate case 

the magnitude of t and i-| is the same, meaning that 

the splitting of the bottom two states with different val- 
ley composition is the same as the splitting at the usual 

anticrossing given by i (i.e. the anticrossing in Figsj2] 

and|3]) 

The analysis presented underscores the importance of 

determining the relative value of t and This can 

be accomplished by setting up the DQD for a resonant 
tunneling experiment, as in Fig. [6] It is necessary first to 
do a transport experiment through each individual dot 
in order to identify the single-dot energy levels. In an 
ideal situation fc^T <C A ^ e^i, and the levels on each 
dot have the form shown schematically in Fig. [6] Sub- 
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Figure 4: Singlet states when A =0.1meV, t =0 and all 

the interdot tunneling is intervalley. On the far right side of 
the graph we recover the one-dot singlets, as in Fig. |2] 
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Figure 5: Triplet states when A =0.1meV, t-- = and all 
the interdot tunneling is intervalley. On the far right side of 
the graph we recover the one-dot triplets, as in Fig. [3] 



sequently, one can do a resonant tunneling experiment 
through the double quantum dot. The first peak in the 
resonant tunneling current will be between the lowest 

two levels and will give t The levels of the left dot 

are subsequently held fixed while the levels of the right 
dot are scanned. The next peak in the resonant tunneling 

current should then give In addition to providing 

values for the tunneling parameters, such an experiment 
should also give a good indication of the length scale 

characterizing interface roughness correlations. If is 

close to zero it indicates that the correlation length of 
interface roughness is either much smaller than the size 
of a quantum dot or much larger than it. On the other 

hand, if is measurable and of comparable magnitude 

to t it indicates that the phase of A is different on 

the two dots and, from our discussion, roughness correla- 
tions occur on a scale comparable to that of a quantum 
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Figure 6: Resonant tunneling in valley split quantum dots. 
The matrix elements for tunneling between like valley eigen- 
states and between opposite valley eigenstates are expected 
to be very different in magnitude. Resonant tunneling can de- 
tect these tunneling rates and give an accurate indication of 
their relative strength, allowing one to infer the length scale 
of interface roughness correlations. 

dot. For this experiment, as for the valley manipula- 
tion scheme proposed in this paper, it is necessary to 
know the value of the valley splitting so that valley-split 
levels can be unambiguously identified. Two potential 
risks are associated with such a resonant tunneling ex- 
periment. Firstly, if the valley splitting is smaller than 
the typical broadening of the current peaks, then the sig- 
nals due to t and t-| may merge into one peak. In 

this case, the current peaks can be expected to have the 
same magnitude between any pair of levels on the left 
and right dots. Secondly, it may happen that one sample 

will have \t | = \t^ | within experimental error. Aside 

from the fact that such a risk is small, it can be over- 
come by performing the experiment on several samples. 
Furthermore, since 2A is expected to be smaller than the 
confinement energy, the two tunneling parameters can be 
unambiguously identified. 

VII. DISCUSSION 
A. Overview of length scale considerations 

Our analysis makes evident the fact that the length 
scale characterizing interface roughness correlations plays 
an important role in determining the effect of roughness 
on the valley-orbit coupling and intervalley dynamics of 
double quantum dots. The two possible scenarios, inter- 
face correlations of the order of the lattice constant and 
greatly exceeding the lattice constant, produce vastly dif- 
ferent physics on a qualitative level, as well as quantita- 
tively different results. 

When the length scale of interface roughness correla- 
tions is of the order of a few lattice constants the mag- 
nitude of the valley-orbit coupling differs slightly on the 
two dots due to the local roughness, yet its phase is the 
same. All the interdot tunneling occurs between like val- 



ley eigenstates, that is, -I- to -I- and — to — , while inter- 
dot intervalley tunneling is suppressed. We distinguish 
between two limiting physical situations depending on 
the strength of the local roughness relative to the inter- 
dot tunneling matrix element. Weak local roughness is 
characterized hy to ^ v for le2d and t ^ v for 2e2d. 
In the le2d setup the lowest energy levels do not cross, 
they come to within a distance to of each other. Since 
to 3> V, the local roughness cannot mix these levels. 

A feasible valley manipulation scheme using the sin- 
glet states cannot be devised either. Firstly, the sin- 
glets are only separated clearly from other levels in the 
(0,2) configuration. Secondly, there is no crossing of sin- 
glet levels that could be exploited. At the points at 
which singlets cross other triplet states are also close by. 
Even though roughness does not mix singlets and triplets, 
other mechanisms, such as the hyperfine coupling to the 
nuclei, could. On the other hand, the lowest 2e2d triplets 
do cross, providing a scheme that is experimentally ac- 
cessible. By applying a magnetic field experiment can 
unambiguously access the T"*" triplet corresponding to 
the setup in which both electrons are on one dot. Once 
this state is initialized the detuning between the dots is 
ramped up until this state crosses the next-lowest energy 
state. At this point the local roughness can mediate tran- 
sitions between these two levels, and sets the time scale 
for this mixing, which will be estimated below. 

For strong local roughness, v ^ i, it is obvious that 
the lowest energy levels cross in all setups. There is a 
tunneling matrix element between the two lowest energy 
levels, which in all three setups have different valley com- 
position. The time scale for this level mixing is set by the 
interdot tunneling, in contrast to the weak local rough- 
ness case, in which it was set by the local roughness pa- 
rameter V. On the one hand tunneling can be adjusted 
by means of the gate that controls the interdot barrier, 
allowing some experimental control. 

When the length scale of interface roughness correla- 
tions is of the order of the size of a quantum dot, the mag- 
nitude of the valley-orbit coupling will again be slightly 
different between the two dots, yet its phase can in princi- 
ple also be different. The relative phase of the valley-orbit 
coupling is determined by the value of the function e^*''°'» 
on each dot and, where 2ko sets a length scale of « 0.5A. 
The relative phase of the valley orbit coupling determines 
the relative magnitude of interdot tunneling between like 
valley eigenstates and interdot tunneling between oppo- 
site valley eigenstates. Given the random nature of C, the 
extremely short length scale set by 2fco and the impos- 
sibility of achieving control over such short distances in 
the foreseeable future, the relative strengths of tunneling 
between like and opposite valley eigenstates under these 
circumstances is in effect a random number, expected to 
vary greatly from one sample to another. 
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B. Strong intervalley coupling limit 

The extreme cases considered, t_| = t and t = 

are rather similar quahtatively. Once again one must con- 
sider how the strength of the local roughness compares to 
the interdot tunneling. For t ^ v the lowest energy levels 
do not cross in any of the three possible setups. If w » t 
the bare Hamiltonian is the same as when the tunneling is 
predominantly intravalley. We have established that the 
lowest energy levels do cross in all setups, yet the tunnel- 
ing between them is determined by the relative strength 
of t and Since this is a random number any val- 
ley manipulation scheme under these circumstances must 
be regarded as overly ambitious. If interface roughness 
correlations occur on the scale of a quantum dot the val- 
leys are mixed randomly and any attempt at coherent 
manipulation of valley eigenstates would be futile. 

One important result of our study is the fact that vari- 
ations in the position of the interface (given by (C^)) and 
variations in the local roughness (given by v) influence 
the valley-orbit coupling and interdot dynamics in vastly 
different ways. The length scale of the interface misalign- 
ment across a double quantum dot gives rise to two quali- 
tative scenarios: short correlations of interface roughness, 
over the scale of a few lattice constants, and long correla- 
tions of interface roughness, over the scale of a quantum 
dot. For short correlations the valley-orbit coupling is 

suppressed by a factor e 2 ^ with k determined by 
the locations of the valleys in reciprocal space, and (C^) 
the variance of the fluctuations in the interface location. 
In the absence of local roughness, the valley-orbit cou- 
pling in a double quantum dot is the same on each dot. 
For long correlations the phase of the valley-orbit cou- 
pling can also be different between the two dots, such 
that interdot tunneling between opposite valley eigen- 
states is enabled. In this case the relative phase of the 
valley-orbit coupling cannot be controlled and the rel- 
ative magnitude of interdot tunneling between like and 
opposite valley eigenstates is essentially a random vari- 
able. These findings are general and applicable to all 
systems where valley physics is important. 



C. Ideal scenario 

The analysis presented in this work suggests that the 
ideal scenario for inducing transitions between valley 
eigenstates is to use the polarized triplet T"*" states in 
a DQD in which t ^ v. To determine the time scale for 
these transitions we need to estimate the magnitude of v, 
or more precisely its variation across a DQD (which may 
be due to dangling bonds or impurities.) In the absence 
of experimental data, it is necessary to estimate the vari- 
ation in Ad between the two dots. We assume =0.1 
meV and, to cover a wider range of experimental possibil- 
ities, we consider variations in A^j of 1 — 10%. As shown 
above, knowledge of the difference in Ad enables one to 



determine the difference in Vzz 
sume wf^ ^ Imv^g in Eq. 24 



between the dots (we as- 
so that the difference in 
Ajj is given by the difference in v^z only.) Moreover, we 
can use the evaluations of Ad and Ad to infer a series of 
facts about the matrix elements of V{x,y). The matrix 
elements Vzz and Vzz follow the same pattern as Ad and 
Ad respectively, except that instead of Uq we have the 
average of V{x, y) over a quantum dot. Therefore, based 
on the ratio of Ad and Ad, we conclude that the ratio of 
the off-diagonal matrix element Vzz to the diagonal ma- 
trix element Vzz can range between « 0.1 — 0.3 depending 
on whether SiGe is used or Si02. We estimate in addition 
the real part of Vzz to be one order of magnitude smaller 
than the imaginary part. The roughness average reduces 
this by another order of magnitude, with a further re- 
duction by an order of magnitude arising from the mul- 
tiplicative factor t/A]j. One ultimately expects v in the 
range « 0.1 — 1 neV, giving a time scale of approximately 
1-10 /iS. This figure represents the time scale for transi- 
tions from one valley eigenstate to another. Gating on 
such a time scale is easily accomplished experimentally.'^ 



D. Important regions of the energy spectrum 

In all scenarios discussed in this work there are three 
relevant regions of the energy level spectrum as a func- 
tion of interdot detuning 6. These are the two extreme 
cases (5 ^ 0, J ^ 0, as well as a third value 6 at which 
the lowest two energy levels cross, in the cases when such 
a crossing does occur. For le2d in the strong tunneling 
regime the two lowest energy levels as a function of de- 
tuning go through an avoided crossing where their mini- 
mum separation is of order t. The matrix element mixing 
these levels is of order v, and since v <^ t the possibility 
of mixing these levels by means of interface roughness 
is remote. For 2e2d in the strong tunneling regime, the 
lowest singlet states also anticross and cannot be mixed, 
yet two triplet states with different valley compositions 
cross each other and can be mixed by roughness. In fact 
transitions between triplet states offer the best prospect 
of experimental manipulation of the valley degree of free- 
dom using interface roughness. For strong roughness all 
the above scenarios {le2d, 2e2d singlet and triplet states) 
are very similar. Without interdot tunneling many en- 
ergy levels cross, and in the neighborhood of the crossing 
point the lowest two energy levels can be mixed by t. 
Given that the tunneling parameter can be adjusted by 
means of the gate that controls the barrier between the 
dots, it can in principle be reduced to a time scale slow 
enough that a coherent valley experiment is possible with 
currently achievable gating times. 



VIII. SUMMARY 

Interface roughness has two principal effects on quan- 
tum dots made in materials with multivalley energy spec- 
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tra. Firstly, it leads to a suppression of the valley-orbit 
coupling by a factor e~^'°o(C >^ where (C^) is the vari- 
ance of the fluctuations in the position of the interface. 
Secondly, in a double quantum dot it enables transitions 
between valley eigenstates. Such transitions occur in the 
dynamics of one and two electrons in a DQD. In the ma- 
jority of processes discussed, one electron initialized into 
a particular valley eigenstate on one dot tunnels into the 
opposite valley eigenstate on the adjacent dot when the 
detuning between the dots is swept to a certain point of 
the energy spectrum and appropriate mixing time is al- 
lowed. The position of the mixing point and the value of 
the mixing time must be determined individually for each 
configuration - leld, 2e2d singlets and 2e2d triplets, and 
depends also on the relative magnitude of the local inter- 
face roughness and the relevant interdot tunneling matrix 
element. Based on our study of all possible configurations 
we conclude that the likeliest system to observe this phe- 
nomenon is provided by the triplet states in the 2e2d 
configuration, which may be reliably initialized and are 
estimated to mix on a time scale of 1-10 /xs. It must 
be pointed out that, although the possibility of induc- 
ing transitions between valley eigenstates by means of 
interface roughness discussed in this article is one way 



of manipulating valleys, the scheme we have presented 
is not tunable. Instead it must be regarded as a first 
step in the attempt to achieve coherent control of valley 
eigenstates. 
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Appendix A: Valley-orbit coupling 

The main contribution to An comes from the interface 
potential. For a perfectly smooth and perfectly sharp 
interface this is 



/>oc poo poo 
/ dxdy\<j>{x,y)\^ dzQ{-z)\i;{z)fe-^"'''ulir)u,{r) 
-oo^— OO J —CO 

/ dxdy\^{x,y)\^e''^^-^ dze'^'\ 

j^l^ -oo J —OO J —CO 

where = ki, + iQz — 2ik(^. We can approximate the x and y integrals by (5-functions of Q^, Qy, yielding 



(*,|V|*,)« UoN^z^ (A2) 



The matrix element between the same valley eigenstates is given by 

/oo poo poo 
/ dxdy\(P{x,y)f dzei-z)mz)\^\u^ir)f 
-OO J —oo J —oo 



/oo pOO P\J 
/ dxdy\ct>{x,y)fe''^--^ dze^'^-+'' 
J\.%^ -OO ^— CXD J —OO 

In the presence of an interface ip{z) — )■ 'ijj{z — C)- presence of interface roughness 

/oo pOO pOO 
/ dxdy[Uo + V{x,y)]\cf>{x,y)\' / e[-(^ - C)] - 01' e"'"'^^ (r) 
-CO J —CO J ~oo 

/CO /* CO /* oo 

/ dxdy[Uo + V{x,y)]\^{x,y)\^e''i--^ e[-(^ - C)] - 01 

-OO''— OO J — oo 



(Al) 



(A3) 



2 gi(Q,-2fe,)z_ 



(A4) 

The z-integral is performed by first changing the variable of integration to z' = z — (let also i^^ = Qz — 

/OO poo pO 

dzQ[-{z-C)]\ij{z-C)\'^e'''i' = / dz'e{-z')\ij{z')fe'''(^''+^^ ^ dz' \i}{z')\'' e'^^'' . (A5) 

-OO J — oo J — oo 
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We always have z' <0 and tp{z') = N zq e""^ , so the integral is trivial 



/ dz'\^{z' 

J — oo 



(A6) 



Appendix B: Matrix elements of roughness in the 
singlet and triplet manifolds 



The singlet eigenstates are 



We will first give the matrix elements of the rough- 
ness potential between the bare singlet and triplet states. 

These are in fact the same, so {S\V\S) = {f\V\f) for all 
the cases below. There are no matrix elements of the 
roughness potential connecting singlet and triplet states. 
The random potential is written as V{ri) + V{r2)- The 
intrabranch matrix elements are 













{si'i\v\si'i) = 




















0. 





The interbranch matrix elements are 



(Bl) 



(B2) 



?<2 



+ 2*2 



1 

71 



+ S. 



RR 
±± 



/2 V 4 



LR 



+ s. 



RR 
±± 



* ^LR . i oLR 



t X 



LR 



^RR 



(C2) 



Two triplet eigenstates are T^±. The other three are 



All other matrix elements are assumed to be negligible. 
Appendix C: Singlet and triplet eigenstates of 2e2d 
for V 



We will use the following notation for the energy eigen- 
states 



4 



4 



-6 - VS"^ + 8f2 



±2Ao 



±2Ao 



~6 + VP + 8*2 



(CI) 



T>_ 



T<_ 



rpsym 



4' + 2*2 



^<2 



2*2 



* rf^LR 

^0 



* ^ 



rpLR _i rriti 

7> + 



RR 



'0 



^0 ^0 



(C3) 



-5 - + 8*2 
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